Universal magnetic and structural behaviors in the iron arsenides 
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Commonalities among the order parameters of the ubiquitous antiferromagnetism present in the 
parent compounds of the iron arsenide high temperature superconductors are explored. Addition- 
ally, comparison is made between the well established two-dimensional Heisenberg-Ising magnet, 
K 2 NiF 4 and iron arsenide systems residing at a critical point whose structural and magnetic phase 
transitions coincide. In particular, analysis is presented regarding two distinct classes of phase tran- 
sition behavior reflected in the development of antiferromagnetic and structural order in the three 
main classes of iron arsenide superconductors. Two distinct universality classes are mirrored in 
their magnetic phase transitions which empirically are determined by the proximity of the coupled 
structural and magnetic phase transitions in these materials. 
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I. INTRODUCTION 

Understanding the fundamental properties of the un- 
derlying magnetic and structural order in the newly 
discovered iron-based high temperature superconductors 
(high-T c ) is an essential step in the eventual resolution 
of magnetism's role in the superconducting pairing in 
these systems. The key unifying feature between the iron 
arsenide class of high-T c superconductors and the well 
known cuprate high-T c s is the universal presence of an 
antiferromagnetically (AF) ordered state in close prox- 
imity to the development of superconductivity within 
their respective phase diagrams 1-3 . The undoped, parent 
systems of both classes of high-T c s exhibit long-range, 
AF order 1 ' 4 that is suppressed upon doping and either 
vanishes 5 or weakly competes 6,7 with the onset of super- 
conductivity. While previous exploration of critical prop- 
erties intrinsic to the static AF order in the cuprates led 
to fundamental insights regarding the interactions gov- 
erning the spin behavior in those systems 8 , experiments 
probing the detailed behavior of magnetism in the iron 
pnictide systems are only just beginning. 

The ordered spin structures of the LaFcAsO(llll) 4 , 
NaFeAs(lll) 9 , and BaFe 2 As 2 (122)-type 10 iron arsenide 
compounds possess a common unidirectional antiferro- 
magnetic ordering as shown in Fig. 1 where the AF prop- 
agation vector points both along the long-axis within the 
basal plane and along the out-of-plane, c-axis. The in- 
plane component of the AF wave vector is determined by 
an accompanying or preceding structural distortion from 
tetragonal to orthorhombic symmetry in which the a-axis 
within the iron layers is slightly elongated. While the 
exchange couplings and spin dynamics intrinsic to this 
magnetic structure in the new class of iron-based high- 
T c s have been a topic of considerable focus, relatively lit- 



tle attention has been given to any precise investigation 
of the magnetic phase transition behavior in these sys- 
tems. Indeed, within the bilayer AEFe2As2 class of pnic- 
tides (where single-crystals are most easily synthesized) a 
number of early studies suggested first-order behavior in 
the onset of the AF phase 11-14 . While these reports seem- 
ingly precluded any subsequent study of critical behavior 
in these systems, separate reports also showed continu- 
ous, second-order magnetic phase transitions in a number 
of the same materials 15-17 . Recently, new single-crystal 
studies have appeared showing continuous phase transi- 
tions in a number of iron pnictide variants 7,18-20 thereby 
suggesting that improvements in sample quality and di- 
versity may now facilitate a reliable comparison between 
systems. 

In this article, we present our analysis of the mag- 
netic and structural order parameters in the three main 
classes of iron pnictide superconductors: AFe 2 As 2 (A- 
sitc = Ba, Sr, Eu, ...), AMFeAs (AM = alkali metal = 
Na, Li), and REFeAsO (RE = rare earth = La, Pr, Ce, 
...). Undoped systems within the bilayer, 122 class ex- 
hibit a two-dimensional (2D) Ising-like magnetism whose 
order parameters are shown to mirror closely the spin 
behavior in the prototypical 2D-Hcisenberg-Ising system 
K 2 NiF 4 . Once the 122 systems are perturbed and the 
structural and magnetic phase transitions no longer coin- 
cide in temperature, a fundamentally different magnetic 
behavior emerges. In this regime, when Tjy <Ts, a dra- 
matically altered development of the magnetic order pa- 
rameter appears and collapses onto a seemingly universal 
curve with a critical exponent of j3 « 0.25. This reflects 
a crossover from the coincident T^ =Ts, 2D-Ising-likc 
behavior to a distinctly separate universality class where 
the coupling between Tg and Tjv has been substantially 
weakened. Whereas the previously observed coupling be- 
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FIG. 1: Crystal structures of representative materials from 
the (a) REFeAsO, (b) AMFeAs, and (c) AFe 2 As 2 classes of 
iron arsenides. Boxes show the enclosed chemical unit cells of 
each system. To the right of each chemical structure are plot- 
ted the ordered spin systems. Moment directions are plotted 
as yellow arrows located on the Fe-sites. 



tween the primary structural and magnetic order param- 
eters in the undoped 122 systems renders identical criti- 
cal exponents for both phase transitions, surprisingly, we 
find that far from the critical point (where Tjy <Ts) both 
critical exponents modeling the structural and magnetic 
phase behavior in the iron arsenides change in an iden- 
tical fashion. Thus even for systems such as LaFeAsO, 
where T s =l.lTjv, the structural, orthorhombic distor- 
tion develops in a manner identical to the onset of AF 
order in the material. These empirical observations of 
potentially universal magnetic and structural behavior 
among the various classes of iron pnictide high-T c s hold 
promising implications for the identification of a common 
role for magnetism within their superconducting proper- 
ties. 



II. MAGNETISM IN K 2 NiF 4 AND BaFe 2 As 2 

Our previous investigation of AF order in the 
BaFe2As2 (Ba-122) compound 18 revealed that the mag- 
netic order parameter was well modeled by a simple 
power law behavior of (1 — with an exponent of 

(3 = 0.125 thereby suggesting that the universality class 
of the AF phase transition was the same as that of the 
2D-Ising model. While a weakly first order component to 
the phase transition could not be entirely precluded, our 
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FIG. 2: Comparison plots of the magnetic order parameters 
in BaFe 2 As 2 and K 2 NiF4 with normalization described in the 
text. The solid yellow line denotes the expected behavior 
for the ideal 2D-Ising system, (a) Neutron measurements of 
M(^) for BaFe 2 As 2 and K 2 NiF 4 taken from Refs. 18 and 21 

respectively, (b) Plot showing M 2 (^-) for BaFe 2 As 2 fit with 
a power law model convolved with a Gaussian distribution of 
Tats as described in the text. 



previous analysis placed a limit on any possible linear, 
first order trade-off between phases to less than 0.5K of 
Tjv. The remarkable range in reduced temperature over 
which the simple power law models the magnetic phase 
transition in BaFe2As2 is reminiscent of phase behavior 
in known low dimensional magnets such as the prototyp- 
ical 2D magnet K 2 NiF 4 21 which stands as a useful stan- 
dard for comparison when examining the phase transition 
behavior in BaFe2As2- 

K2N1F4 is a model 2D system with a crossover from 2D 
Heisenberg behavior at low temperatures to anisotropy- 
driven 2D Ising behavior in the vicinity of Tjy. Figure 
2 (a) displays a comparison of the magnetic order pa- 
rameters, M, as a function of (1 — T/T/v) for both sys- 
tems K 2 NiF 4 and BaFe 2 As2 where data for BaFe2As2 
are taken from Ref. 18 and the data for K2MF4 are 
extracted from Ref. 21. In order to facilitate a compar- 
ison of the phase transition behavior in each respective 
system, all data were normalized at (1 — (T/T/v)) = 0.1, 
where w 75% of the saturated moment had been reached. 
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Examination of this comparison highlights the striking 
similarity between the evolution of the magnetic order 
parameters in Ba-122 and K 2 NiF 4 where the critical be- 
haviors of both systems with freely refined exponents of 
j3 = 0.10 and j3 = 0.14 respectively bracket the ideal 2D 
Ising curve with f3 — 0.125. In particular, this demon- 
strates explicitly that the sharp onset of antiferromag- 
netism in BaFe2As2 mirrors that of a known second or- 
der 2D-Heisenberg-Ising AF phase transition and thus 
that it cannot be easily dismissed as an artifact from an 
obscured first order phase transition. 

At temperatures immediately above Tjv, appreciably 
more magnetic scattering at the 3D ordering wave vec- 
tor persists in BaFe2As2 relative to K2MF4. The ori- 
gin of this tail of scattering above Tjv is likely due to a 
small distribution of AF ordering temperatures through- 
out the bulk of the sample. This can be confirmed by 
fitting the expected power law behavior of the order 
parameter weighted by a Gaussian distribution of T/v's 
within the sample. The result of such a fit to the form 

M(T) 2 = B 2 f (1 - ^-)2/3 e -(T-(T„»/2 CT 2 ig ghown in Fig 

2 (b), and the fit gives values of (T) = 136.12 ± 0.03K, 
(3 = 0.106 ± 0.002, and a = 0.28 ± 0.02K. The power 
law was forced to for T > TV, and only data in 

the range 0.85K< < 1.025K were included. The 

N 

excellent agreement between this Gaussian-broadened 
power law description and the magnetic order param- 
eter in BaFe2As2 clearly demonstrates that the phase 
transition in its entirety can be effectively modeled by 
a simple power law with a (3 w 0.125 and a slightly 
broadened T/v with full width at half-maximum, FWHM 
= 0.0048 in reduced temperature. Additionally, our sub- 
sequent fitting of the same data to the first order form 
M 2 (T) = (A s +B(l- ^Y) 2 yielded negligible A s values 
and identical (3 values to the previous A s = fits, thus 
further supporting our claim of continuous behavior in 
BaFe 2 As 2 . 



III. MAGNETIC AND STRUCTURAL ORDER 
PARAMETERS IN THE IRON ARSENIDES 

Having established that the magnetic phase transi- 
tion in Ba-122 mirrors that of a known 2D-Heisenberg- 
Ising magnet, it is informative now to compare the mag- 
netic order parameters between different classes of iron 
arsenides. In Fig. 3 (a), measurements of magnetic 
order parameters squared are plotted for both Ba-122 
and LaFeAsO, where an immediate difference between 
the temperature dependencies of the order parameters is 
apparent. The data for LaFeAsO are taken from Ref. 
22 which also reports a clear commonality among the 
measured magnetic order parameters of the REFeAsO 
systems. Therefore from Fig. 3 (a) it is clear that 
there is a fundamental difference in the phase transi- 
tion behavior of the Ba-122 and the 1111 class of iron 
arsenides where the order parameter in LaFeAsO fit be- 
tween 0.6 < -Sr- < 1.05 to a simple power law gives a 
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FIG 

BaFe2As2 and LaFeAsO. The solid black line denotes the 
Gaussian convolved power law fit from Fig. 2 (b), and the 
dashed line shows a power law fit with j3 = 0.25 as described 
in the text, (b) Overplot of the magnetic order parameter 
(data taken from Ref. 22) and the structural order parameter 
(data taken from Ref. 23) for the LaFeAsO system. 



(3 = 0.22 ± 0.06. The dashed line in Fig. 3 (a) shows 
a simple power law function with a fixed j3 = 0.25 over- 
plotted with the LaFeAsO data. Unfortunately, due to 
the lack of detailed data close to the phase transition in 
LaFeAsO, the critical exponent for the magnetic order 
parameter cannot be determined precisely; however from 
Fig. 3 (a) it is apparent that a simple (3 — 0.25 provides 
a reasonable approximation to the magnetic phase be- 
havior in this system. This renders a critical exponent (3 
in the monolayer iron arsenides roughly two times larger 
than the (3 = 0.125 reported to model the phase behavior 
in Ba-122. 18 Given that both systems possess the same 
intrinsic spin structures and magnetic ions, this sharp 
distinction is surprising and, we believe, of fundamental 
importance. 

In noting that the strongly coupled structural and 
magnetic order parameters in Ba-122 are known to show 
an identical temperature evolution, 18 it is clearly impor- 
tant to examine how the structural phase transition is 
modified in LaFeAsO where T5 and Tn are largely de- 
coupled. An overlay of both the magnetic order param- 
eter data again taken from Ref. 22 and the structural 
order parameter data extracted from Ref. 23 is plotted 
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in Fig. 3 (b). From this figure, it is immediately evi- 
dent that, despite the splitting of Tg and by over 
more than 14K, the phase behaviors of both the mag- 
netic and structural order parameters in LaFeAsO track 
one another quite precisely. The yellow solid line in Fig. 
3(b) shows the results of a power law fit with j3 — 0.25 
overlaid with both sets of data while the black line shows 
the results from the freely refined fit of the magnetic or- 
der parameter in LaFeAsO (discussed previously). This 
simple, empirical, comparison suggests that the magne- 
tostructural coupling observed in Ba-122 and other un- 
doped, bilayer pnictides is also reflected in the monolayer 
pnictides where Ts and Tjy are no longer coincident. The 
primary effect of tuning away from the Tjv=Ts critical 
point appears to be a renormalization of the critical ex- 
ponents modeling both the structural and magnetic order 
parameters which are altered in the same fashion. De- 
spite the significant difference in (3 values between the 
122 and 1111 compounds, the structural and magnetic 
phase behaviors mirror one another identically in both 
systems, thus strongly suggesting that this is a universal 
feature within the parent phases of the iron arsenides. 

In considering the immediate contrasts in the mag- 
netic phases of the 122 and 1111 compounds which may 
lead to the large differences in their magnetic and struc- 
tural phase transition behaviors, the primary distinction 
is that the strongly coupled, concomitant, magnetostruc- 
tural phase transitions in the 122 series splits in tem- 
perature within the 1111 compounds thereby potentially 
altering the critical spin behavior. In order to investi- 
gate this further, we examined data from the literature 
reporting the magnetic order parameters in lightly Ni- 
and Co-doped Ba-122 systems 7,19 ' 24 . Within these ini- 
tial reports, it was observed that doping small amounts 
of electrons or holes splits the magnetic and structural 
phase transitions in Ba-122, and from the reported data 
it can be clearly seen that the observed order param- 
eter is dramatically modified from the parent system's 
behavior. In both cases, a simple power law with an ex- 
ponent of (3 ~ 0.3 was reported to model the magnetic 
phase transitions close to TV 7,19 - In order to provide a 
more explicit comparison between these various systems, 
Fig. 4 overplots the magnetic order parameters of both 
of these doped Ba-122 systems along with the magnetic 
order parameters of several different iron arsenide parent 
systems extracted from the literature 7 ' 9 ' 18,19 ' 24,25 . Sys- 
tems whose phase transitions are closely modeled by a 
j3 = 0.125 were normalized together at T/Tn = 0.9 fts 
0.56M 2 (^- = 0), and systems whose phase transitions 
were approximated by (3 ~ 0.25 were normalized together 
at T/T N = 0.68 w 0.56 M 2 (^ = 0). For clarity, the or- 
der parameters of each class of systems were set apart 
by cross-normalizing the average values of the order pa- 
rameters of the (3 « 0.25 class and the (3 — 0.125 class 
at 75% of their saturated values. In determining which 
data to include in this comparison, we attempted to se- 
lect the highest quality data in the literature using single 
crystal measurements when possible. As Sn-flux grown 




FIG. 4: Plot showing the comparison of M 2 (^-) for systems 
in the 122, 1111, and 111 classes of iron arsenide supercon- 
ductors with the normalization procedure between different 
systems described in the text. The solid black line and dashed 
black line denote the ideal 2D-Ising and 3D-Ising values for 
M 2 (^) respectively. The collapse onto two distinct curves 
for systems where Tjv=Ts and for systems where Tjv <Ts is 
shown. 



crystals exhibit Sn incorporation which strongly renor- 
malizes magnetic properties in certain 122 pnictides 26 , 
we utilized data only from studies performed on crystals 
grown via FeAs or self-flux methods. 

Looking at Fig. 4, it is immediately clear that there are 
two distinct classes of transitions amongst the data plot- 
ted. Both 122 parent systems, SrFe2As2 and BaFe2As2, 
collapse onto a single curve that is modeled well by the 
/3 = 0.125 2D-Ising order parameter. Data from doped 
122 systems and from 1111 and 111 iron arsenide com- 
pounds in Fig. 4 reveal that, in systems where the struc- 
tural and magnetic phase transitions separate in temper- 
ature, the resulting magnetic phase transition collapses 
onto a separate universal curve whose critical exponent 
is roughly double that of the Tjv =Ts materials. Within 
this second class of magnetic phase transitions, fits to the 
data between 0.6 < ^ < 1.1 led to an average exponent 
of j3 = 0.27 ± 0.03 demonstrating that, once the struc- 
tural phase transition decouples from the onset of mag- 
netic order, the universality class of the magnetic phase 
transition is clearly altered. Currently however, the qual- 
ity of data from these systems with broader, decoupled 
magnetic phase transitions precludes the determination 
of their precise critical exponents; however we have plot- 
ted in Fig. 4 one possibility of a crossover to a 3D-Ising 
order parameter with (3 = 0.325. A 3D-Heisenberg tran- 
sition with f3 — 0.365 (not shown) may also fall within 
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the error of the data suggesting in both cases an increased 
dimensionality to the magnetic phase transition once de- 
coupled from Tg. 

IV. DISCUSSION 

Several theories have been put forth examining the 
relationship between the higher temperature structural 
transition and the lower temperature magnetic phase 
transition in the 1111 scries of pnicides. Work by Xu 
ct al. considered the presence of a finite temperature 
2D Ising order which can be envisioned as an electronic, 
nematic phase that forms prior to the onset of the spin 
density wave (SDW) order. 27 Additionally recent work 
by Chen et al. has considered the presence of an or- 
bital ordering parameter also possessing 2D-Ising symme- 
try stemming from a simple spin-orbital Hamiltonian. 28 
SDW order in these models freezes out due to the finite J c 
coupling and is simply a reflection of the eventual onset of 
three-dimensional exchange within the system. However, 
to the best of our knowledge, no theoretical work explor- 
ing the expected spin behavior at the effective critical 
point in the 122 pnictides where both phase transitions 
are coincident in temperature has been performed. 

Experimental work on BaFe2As2 has reported a re- 
duction in the anisotropy gap of the spin excitations 
upon Ni-doping onto the Fe-site, potentially indicating 
an enhanced two dimensionality induced upon carrier 
doping. 19 This seems to run counter to our current ob- 
servation of a transition from a regime in which the mag- 
netic phase transition behavior is consistent with that of 
a two dimensional magnetic system in the BaFe2As2 sys- 
tem to a regime in doped samples whose larger critical 
exponents (ie. [3 w 0.27) suggest a new universality class 
with a more three dimensional character. Recent pho- 
tocmission experiments however have shown an abrupt 
change in the dimensionality of Co-doped Ba-122 where 
the two dimensional electronic states of the parent Ba- 
122 system are tuned toward a more three dimensional 
dispersion as Co is doped into the system. 29 Separate 
neutron studies of BaFe 2 As 2 have also reported largely 
two dimensional spin fluctuations above T5 17 consistent 
with the picture of a two dimensional magnetic order pa- 
rameter close to Tn and qualitatively in agreement with 
the largely two dimensional Fermi surface reported in 
recent photoemission measurements on BaFe2As2 29 and 
CaFe2As2. 30 To-date however, there has been no com- 
prehensive experimental investigation into the detailed 
evolution of the dimensionality of spin fluctuations in a 
122 system as it is tuned from the quasi- 2D high tem- 
perature regime above T5 to the anisotropic, three di- 



mensional magnetic and electronic phase far below Tjv- 
Future work exploring this along with further theoreti- 
cal efforts detailing the coupled phase transitions in the 
pnictides are needed in order to understand the seem- 
ingly universal structural and magnetic phase behaviors 
within these materials. 

V. CONCLUSION 

In summary, we have presented an empirical analysis of 
the magnetic order parameters in the primary classes of 
the iron arsenide superconductors. Within the undoped, 
bilayer 122 class of iron pnictides where Ttv^Ts, both 
the magnetic and structural phase transitions are well 
modeled by a 2D-Ising order parameter; however, upon 
doping, the magnetic and structural phase transitions al- 
ter to exhibit a critical (3 ss 0.27. The likely cause of 
this crossover in the phase behavior is the decoupling of 
the magnetic and structural transitions where alternate 
classes of iron pnictides with T^v <Ts collapse onto the 
same seemingly universal curve. Two distinct behaviors 
therefore appear in the magnetic and structural phase 
transitions of the iron arsenides: The first of these oc- 
curs when the structural and magnetic phase transitions 
coincide at a multicritical point, and the second emerges 
once these two transitions no longer coincide rendering 
strongly renormalized critical exponents for both. The 
critical behavior of the iron pnictides therefore seems to 
transition from one which parallels a known 2D mag- 
net when Tjv=T,s to one consistent with a more three 
dimensional character upon tuning away from the criti- 
cal point. This observation currently contradicts several 
existing arguments for more two dimensional behavior 
reported in the Ba-122 system upon introducing charge 
carriers, thus suggesting that new theoretical insight is 
needed. Our analysis demonstrates that further theoret- 
ical and experimental work exploring the phase behavior 
in the critical magnetism of the iron pnictides have the 
possibility of resolving exciting new physics relevant to 
the fundamental symmetries in the magnetism of these 
compounds. 
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